Planck scale inflationary spectra from quantum gravity 
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We derive the semiclassical evolution of massless minimally coupled scalar matter in the de Sitter 
space-time from the Born-Oppenheimer reduction of the Wheeler- De Witt equation. We show that 
the dynamics of trans-Planckian modes can be cast in the form of an effective modified dispersion 
relation and that high energy corrections in the power spectrum of the cosmic microwave background 
radiation produced during inflation remain very small if the initial state is the Bunch-Davies vacuum. 
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I. INTRODUCTION 

It is commonly accepted that inflation lj is a viable 
paradigm for the early Universe which solves some of 
the problems of the standard big-bang scenario and al- 
lows one to make testable predictions about the spectrum 
of the cosmic microwave background radiation (CMBR). 
However, it has also been realized that inflation pro- 
vides a window towards trans-Planckian physics [2| as 
it magnifies all quantum fluctuations and red-shifts orig- 
inally trans-Planckian frequencies down to the range of 
low energy physics currently observed. This has raised 
two related issues: a) how to describe the original quan- 
tum fluctuations in such a high energy regime, given that 
there is currently no universally accepted theory of quan- 
tum gravity, and b) whether the effect of now red-shifted 
trans-Planckian frequencies can be observed with the pre- 
cision of present and future experiments. 

One can take the pragmatic view of the renormal- 
ization group approach to high energy physics and use 
quantum field theory on a classical metric background 
to describe quantum fluctuations after their frequencies 
have been red-shifted below the scale of quantum grav- 
ity. However, the higher energy dynamics cannot be 
a priori neglected as it could enter (at least) in the 
form of initial conditions for the quantum state of the 
field theory 0, Q and, in turn, affect the CMBR spec- 
trum (to an extent which is currently being debated; see, 
e.g., Refs. @SHH0). 

As a consequence of the above observations, one 
expects that the dispersion relations of matter fields 
change on approaching the Planckian regime, so that 
the frequency to of massless modes will depend on their 
wavenumber k according to an expression of the general 
form 



w = - [l + f(k,£ p /a,a,...)] 
a 
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in which £ p is the Planck length and a = a(t) the scale 
factor of the flat Robertson- Walker metric 



ds 2 = -N 2 dt 2 



(dxy 



(2) 



Most attempts have tested the effects of functions / cho- 
sen ad hoc in Eq. (JJJ 2] . We shall instead derive the dis- 
persion relation for a minimally coupled massless scalar 
field from the action principle previously employed to 
study the semiclassical dynamics of non-minimally cou- 
pled scalar fields in Ref. 

In that paper, the principle of time-reparameterization 
invariance was lifted to a quantum symmetry. We 
then obtained an Hamiltonian constraint from which the 
Born-Oppenheimer (BO) reduction ^(j allowed to prop- 
erly and unambiguously recover the semiclassical limit of 
quantum field theory on a curved background starting 
from the Wheeler-DeWitt (WDW) equation 11}. This 
required that certain "quantum fluctuation" terms be 
negligible with respect to the usual matter contributions. 
Precisely such terms in the matter equation will now be 
studied and their effect on the power spectrum derived 
for the simplest model of inflation. 

In the next Section, we briefly introduce the model 
and its classical and Wheeler-DeWitt dynamics in the 
BO decomposition (for more details, see Ref. @). In 
Section [□] we then focus on the quantum equation for 
matter which we solve perturbatively in order to deter- 
mine a dispersion relation and the corresponding CMBR 
spectrum. Finally, we comment on our findings in Sec- 
tion Cvl We shall use units with c = h = 1. 



II. SCALAR FIELD IN DE SITTER SPACE 



As the simplest model of inflation pi 
sider the de Sitter metric with 

a = ciq exp(7i t) 



we shall con- 



(3) 



and cosmological constant A = 3H 2 . The massless min- 
imally coupled scalar field (f> will be treated as a pertur- 
bation with respect to A in the interval t\ < t < if, where 
the relevant modes k were first produced at t = t\ and 
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inflation ends at t — if . Upon varying the corresponding 
action (/ = d t f), 



S = - 



'S ( a 2 
a 3 N dt 1 



1 



21% J u \N 2 a 2 



a A Ndt / d^x 
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~3 
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N 2 
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(4) 



one obtains the Euler-Lagrange equations of motion in 
the form of the Hamiltonian constraint 



S N S = -H = 
and the dynamical equations 

S a S = S^S = . 



(5a) 



(5b) 



The former allows one to work in the gauge N(t) = 1, 
provided the initial conditions are such that H(ti) = 0, 
and Eqs. (|5b|l will then evolve the initial data consis- 
tently, since the Hamiltonian constraint l|5a|) is then pre- 
served, 



B. BO decomposition 

We can decompose the wavefunction of the Universe 
according to the BO prescription [T(| as 



|¥(a,0>=#»)|x(a,0>^]X 



\Xk) 



(11) 



where \ xk) is the wavefunction for the mode k. Eq. © 
can then be shown [l2| to be equivalent to a system of 
decoupled Schwinger-Tomonaga equations (one for each 
mode) and a Friedmann-like equation expressing the 
back-reaction of matter on the expansion of the Uni- 
verse. In fact, starting from the ansatz one can 
reduce Eq. © to the two coupled equations 



f<^)V> (12a) 



-f(D_ij) D + \ X )+^[H^-(H lf> )) \x) 



-^^(d 2 

2a Y V + 
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H = —a5 a S — 4> S^S = . 
However, since 

H{t) =6^ = => 6 a S = 



(G) 



(7) 



for our purposes it is more convenient to impose the 
Hamiltonian constraint together with the Klein-Gordon 
equation at all times. 



A. WDW equation 



At the quantum level, Eq. I|5a(l becomes the WDW 
constraint which reads 



H\*) = (H G + J2 k Hk)m = 0, (8) 
where the gravitational Hamiltonian is 01 



2 V p a 3£* 



(9) 



and the matter Hamiltonian is given by the sum of the 
contributions for each mode <fik of the scalar field, 



H k = U§ + ak 2 4>t, 



(10) 



with [4>k, it k'] = i$k,k'- This equation is rather involved, 
but the following treatment allows us to obtain more 
tractable equations. 



where (O) = (x \0\x) and the action of the "covariant" 
derivatives D ± = d a ± (d a ) on if) and | x) reduces to that 
of d a after rescaling by a geometric phase , 



|x) - \x)^c +ir{d «' )da '\x) 



(13a) 
(13b) 



Eq. I|12b|) is still rather involved since it contains 
the full matter wavefunction. Upon projecting onto 
Iln^fc (X« L one ean obtain independent equations for 
each mode, 



da\Xk) + (Hk-(H k ) k )\ X k) 



\Xk) 



(14a) 



where (0) k = (\k |0| Xfe) and the terms in the right hand 
side represent the kind of "quantum gravitational fluctu- 
ations" already mentioned in the Introduction. Finally, 
one can also write the gravitational Eq. I|12a|) in terms of 
the solutions to the matter Eq. I|14a|l as 



^d 2 a 
I a 



"' A E<^)^ = -^E(^(i4b) 



6 11 



Again, "quantum gravitational fluctuations" appear in 
the right hand side which will affect the way matter back- 
reacts on the evolution of the scale factor. However, we 
shall not consider the latter equation any further in the 
present paper. 
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III. MATTER EQUATION 



A. Perturbative analysis 



Our aim is to solve Eq. (|14af> for modes with wave- 
lengths a/k < l p <C H.~ 1 at the time t ~ t^ and evolve 
them to a later time t = tf when a/k ^> Hr 1 in order to 
determine the power spectrum. The choice of the state 
| vac) is very important and will be discussed later. First 
we take the semiclassical limit for gravity which allows 
us to introduce the time as 0, Q 



( 2 



(d a In ipj d a 



iadt 



(15) 



In the infrared sector, a/k ^> £ p , the right hand side 
(RHS) of Eq. 114a|l can be discarded, since 



RHS ~ fa) d 2 a \ xk) ~ £ (k'/a 3 )] Xk ) 



(16) 



is suppressed by a factor of order k£ 2 /a 2 <C 1 with re- 
spect to (Hk) ~ k/a. Eq. Q14a|l therefore coincides with 
the Schwinger-Tomonaga equation for \xk) 0, El 0|- 
For the same reason, however, one expects that the RHS 
of Eq. (|14af> significantly affects the scalar field dynamics 
in the ultraviolet sector a/k < £ p . If we define 



\ Xs )=e- i ! t ^« dt '\x k ) 



Eq. I|14al) then becomes 
3i5 2 



1 



2a 3 H 3 
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Hk)\Xs 



2a 3 H 3 



(17) 



A|Xs) ,(18) 



where 6 = H£ p and A = £ 3 =1 (d t - ( Xs | 6 { \ Xs )) with 

6 X = 2H- 1 ( Xs \H k \x s )idt (19a) 
d 2 = 3iH k (19b) 
d 3 =H- 1 d 2 . (19c) 

Eq. (|18fl is not linear and the superposition principle no 
longer holds. However, we note that in a (semi) classical 
Universe the dimensionless parameter S 2 <C 1, and we 
can therefore identify two regimes: I) in the very early 
stages a 3 Ti 3 < 8 2 matter evolves according to 

(i d t - H k ) I Xs) - ^ £). A * I Xb> = Wj I x s ) ; (20a) 

II) after a has become sufficiently large (a ^> 8 2 / 3 Tl~ 1 ), 
Eq. l(T%|) can be expanded to leading order in 5 as 



(id t ~ H k ) \ xs 



2a 3 H 3 



J2^\Xs)^S 2 W n \Xs)(20b) 



These expressions represent our main qualitative result: 
Eq. (|20bl) shows that, at late stages of the cosmologi- 
cal evolution, corrections coming from W\\ are of order 
8 2 -C 1 and thus very small; further, although the RHS 
of Eq. (|14a|l seems to produce corrections of order 6° ~ 1 
in the very early stages, we shall see that the effect of W\ 
is actually negligible in Eq. (|2Ua[l . 



On neglecting W (= W\ or W\\), the matter equa- 
tions l|2(Ja|) and l|20b|l simply read 



k ( 1 

id t \Xs) = H k \x s ) = - a f a+- 
a \ I 



where 



k a 2 



(21) 



(22) 



and are the usual annihilation and creation operators 
such that [a, a^] = 1 and can be used to evaluate the 
energy of the matter state. 

Because of the time dependent a = a(t), Hamiltonian 
eigenstates defined by a^ a \ tie) — n \ tie) do not satisfy 
Eq. (|21() . A basis of exact solutions of the time-dependent 
problem is instead given by eigenstates of the invariant 
number operator b I n) = n \ n) , where |20j | 



V2 



H>k i ■< 

Yi (p TTfe - a p4> k 

p v 



and p = pit) must satisfy 



p + 3Hp- 



k 2 



a z 



a 6 p 3 



(23) 



(24) 



One then has [b , = 1 and exact solutions are given by 
superpositions of the base vectors 



|n) 



Ji n G 



-(& + r|o) 



with b 1 0) =0 and the phase 



e(t) = 



at' 



o a 3 (t')p 2 (t') 



(25) 



(26) 



Note that the invariant operators b and w (as well as 
the states |n)) are a mathematical tool to determine the 
solutions and do not in general have a physical meaning. 
They are however related to d and a^ by the Bogoliubov 
transformation [l7| 



= B* b + A* P 



a Vk p + 
a Vk p 



1 



a Vk p 
1 
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(27) 



which also relates Hamiltonian eigenstates to exact solu- 
tions. In terms of such operators, one finds 



W\n) 



12 * 

ab —a 



+f3b 4 +f3* (b^y 



\n) (28a) 



4 



where terms proportional to f3 will be omitted from now 
on since they do not affect the final power spectrum [to 
order S 2 , see Eq. I)44[l] . and 



an 



3 ai 



2 i a 3 n 3 a 5 TL A 



2AB* 



2iH 



\ A \ +o +TTZ 



in 

k/a 



k/c 



(B*y + a 2 



.(28b) 



Since a\/{Hk) ~ aH/k < \S\/k -C f, it appears that W\ 
is actually negligible with respect to Hk- In the regime I, 
one can therefore choose suitable initial conditions p(ti) 
and p(t[) for which 



A(ti) = and B(U) = 1 



(29) 



so that \ue) — \ n) and |xs) = 1"-) &t t = t\. This state 
will not change significantly until it enters the regime II 
and starts to be acted upon by W\\. We can then assume 
a perturbative expansion of the form 



\ Xs ) = \n s )^\n)+8 2 \nW) = M 
From this it then follows that 



i d t R n - 



Hk,Rn 



which admits the solution 



it. In) 



rb z 



>.b 2 - 



S 2 R n 



(30) 



I") ,(31) 



provided r satisfies the differential equation 
if + 28r = a . 



(32a) 



(32b) 



Although R n contains Hermitian parts, matter evolution 
is unitary in general 0| due to the form of Eq. (|18|l and 
in the case at hand one can in fact show that (XslXs) = 1 
to order 5 2 at all times for the states l{3Tj|) . 



B. Dispersion relations 

Two quantum states | Xs) and | Xs) are physically indis- 
tinguishable when they share the same expectation values 
for all the measurable observables X, 



(Xs \X\xs) = (Xs I ^ I Xs) • 



(33) 



In this sense, since we are interested in observable quan- 
tities, such as the power-spectrum, quadratic in <j> and 
7T0, we cannot distinguish between \ n s ) and the state 



1 + iS 2 H r 



where 



H„ = i 



n 2 + n + 1 
2n+ I 



rb — r' 



(34a) 



(34b) 



which evolves in an explicitly unitary way [T^. 

One can for instance consider the invariant vacuum, 
|n s ) = |0) at t = ti, corresponding to an initial state 
devoid of particles (of wavenumber k) . The Hamiltonian 
associated with its effective evolution (|34all is given by 



h: 



eS 



H k -6 2 



ab 2 



+ 7 ( 4>k TTk + TTk 4>k ) ,(35) 



where w is the effective frequency and, for £ 2 /(a 3 H) <C 1 
and aTC/k -C 1, the effective mass is 



/' 



aH 
2a 3 H l~fc 



and 



k 



(36) 



(37) 



Apart from the squeezing term proportional to 7, the 
perturbed dynamics can be described by a modified dis- 
persion relation QJ. In fact, Eq. (|24|l with the initial 
conditions l|29() at t\ — > — 00 admits the solution 



P(t) 



from which one obtains 




321 h 2 e 5 p 

64 a 3 



a 
~ki, 



(38) 



(39) 



valid for trans-Planckian modes with a/k£ p < 1. The 
correction inside the square brackets is proportional to 
the factor (£ p /a) 3 . In a classical Robertson- Walker Uni- 
verse there is no fundamental length scale and one would 
not expect any dependence on the numerical value of 
a = a(i), since only its ratio at two different times has 
a physical meaning. However, in the quantum theory £ p 
plays the role of the fundamental length which breaks 
scale invariance and in terms of which all quantities with 
length dimension must be measured. In fact, "quantum 
gravitational fluctuations" in the right hand side of the 
matter Eq. I|14a|l are precisely proportional to (£ p /a) 3 
and {H£ p ) 2 . 

C. CMBR spectrum 

We can also compute the power spectrum 

7^(fc) = fc 3 (vac|(g|vac) , (40) 

and compare to the well-known case (5 = 0. The re- 
sult will in general depend on the initial quantum state 
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we choose at t = t\ (see, e.g., Ref. As we did for 
the dispersion relation, we choose the adiabatic invariant 
vacuum 



I vac) = 1 S 



(41) 



which initially coincides with the Bunch-Davies vac- 
uum an d is known to give a flat, scale invariant spec- 
trum for modes well outside the horizon, a/k 3> JiT 1 . 
Adding the leading order corrections from Eq. I|28a|) 
yields 

^(fc)~A: 3 ((0 s |<g|0 s )+2<S 2 Re(0 s |$£|0 8 >) ,(42) 
where 



<0 S |^i?|0 s 



(B*) 



A 2 + 2B*A 



(43) 



and r — r(t) is a solution of Eq. I|32b(l . We again employ 
Eq. (|38[l for trans-Planckian modes with k > a(ti)/£ p to 
evaluate the expression 142(1 at t = t{ such that a(t{) ^S> 
k/Ti., and obtain 



r<f,(k)^-H 2 



v \ 



(44) 



where u is an integration constant related to the initial 
conditions for r — r(t) and v ~ 1.22. The zero order 
is precisely the expected flat spectrum, and /c-dependent 
deviations are of order S 2 , as we had already anticipated 
from the general form of Eq. 1201 



IV. CONCLUSIONS 



5 2 = (H£ p ) 2 and thus quite difficult to observe. How- 
ever, in order to obtain the expression in Eq. Q44|). we 
used the condition that the initial state be the Bunch- 
Davies vacuum l@TJ. Different results would follow from 
other choices, such as the one considered in Ref. [3j which 
produces corrections of order TL £ p and thus possibly ob- 
servable 

Si, 

and the non-linearity of the matter equa- 
tion i|18|) may then lead to more non-trivial effects. Our 
result l|44(l can thus be considered as a lower bound for 
the size of quantum gravitational corrections. 

With regard to the two main issues mentioned in the 
Introduction, we can therefore say that our work has 
shown: a) that the semiclassical approach to the WDW 
equation of Refs. can consistently describe the evo- 

lution of trans-Planckian matter modes starting from the 
very early stages and b) that the initial quantum mat- 
ter state can be more relevant than the corrections pro- 
duced during its subsequent semiclassical evolution (see 
also Ref. |j| in this respect). It appears then that ini- 
tial conditions other than 1|41[) could only be justified by 
appealing to new physics beyond the semiclassical level 
of quantum gravity. Quite remarkably, earlier compu- 
tations in the low energy limit of some string and M- 
theories also yielded corrections of order (TC£ P ) 2 for the 
CMBR spectrum |6(, thus in agreement with our (lower 
bound) estimate. 

Let us end by mentioning that our present conclusions 
hold for minimal coupling between matter and gravity 
and may change in more general cases 01 • Moreover, 
the backreaction of the quantum fluctuations in the grav- 
itational equation i|14bjl . which was presently neglected, 
is also worth analyzing. 



We have derived an effective modified dispersion rela- 
tion l|39|) for the trans-Planckian modes of a minimally 
coupled massless scalar field in the de Sitter space-time 
from the action principle and the semiclassical treat- 
ment of the WDW equation 0, ^j. The effects in- 
duced by such modified dynamics on the power spec- 
trum (|44|) during inflation appear as corrections of order 



Acknowledgments 

We would like to thank A.O. Barvinsky, A.Yu. Ka- 
menshchik, G.P. Vacca and G. Venturi for comments and 
suggestions. 



[1] A.R. Liddle and D.H. Lyth, Cosmological inflation and 
large-scale structure (Cambridge University Press, Cam- 
bridge, England, 2000). 

[2] J. Martin and R.H. Brandenberger, Phys. Rev. D 63, 
123501 (2001). 

[3] U.H. Danielsson, Phys. Rev. D 66, 023511 (2002). 

[4] V. Bozza, M. Giovannini and G. Veneziano, JCAP 0305, 
001 (2003). 

[5] A. A. Starobinsky, JETP Lett. 73, 371 (2001); R.H. Bran- 
denberger and J. Martin, Phys. Rev. D 71 (2005) 023504; 
J. Martin and C. Ringeval, JCAP 0608, 009 (2006), 
JCAP 0501, 007 (2005), Phys. Rev. D 69, 083515 (2004); 
H. Collins and R. Holman, Phys. Rev. D 74, 045009 
(2006), Phys. Rev. D 71, 085009 (2005); B. Greene, 



M. Parikh and J.P. van der Schaar, JHEP 0604, 057 
(2006); A. Ashoorioon, J. L. Hovdebo and R. B. Mann, 
Nucl. Phys. B 727, 63 (2005); S. Hossenfelder, Mod. 
Phys. Lett. A 19, 2727 (2004); L. Sriramkumar and 
T. Padmanabhan, Phys. Rev. D 71, 103512 (2005); J. de 
Boer, V. Jejjala and D. Minic, Phys. Rev. D 71, 044013 
(2005); G. Calcagni, Phys. Rev. D 70, 103525 (2004); 
S. Shankaranarayanan and L. Sriramkumar, Phys. Rev. 
D 70, 123520 (2004); R.G. Cai, Phys. Lett. B 593, 
1 (2004); M. Porrati, Phys. Lett. B 596, 306 (2004); 
M. Fukuma, Y. Kono and A. Miwa, Nucl. Phys. B 703, 
293 (2004); K. Ke, Int. J. Mod. Phys. A 20, 4331 (2005); 
S. Hannestad, JCAP 0404, 002 (2004); Q.G Huang and 
M. Li, Nucl. Phys. B 713, 219 (2005); J.P. van der 



6 



Schaar, JHEP 0401, 070 (2004); G.L. Alberghi, K. Gold- 
stein and D.A. Lowe, Phys. Lett. B 578, 247 (2004); 
G.L. Alberghi, D.A. Lowe and M. Trodden, JHEP 9907, 
020 (1999); M. Fukuma, Y. Kono and A. Miwa, Nucl. 
Phys. B 682, 377 (2004); O. Elgaroy and S. Hannestad, 
Phys. Rev. D 68, 123513 (2003); E. Keski-Vakkuri and 
M.S. Sloth, JCAP 0308, 001 (2003); M.B. Einhorn and 
F. Larsen, Phys. Rev. D 68, 064002 (2003); S. Cremonini, 
Phys. Rev. D 68, 063514 (2003); 

[6] N. Kaloper, M. Kleban, A.E. Lawrence and S. Shenker, 
Phys. Rev. D 66, 123510 (2002). 

[7] G.L. Alberghi, R. Casadio and A. Tronconi, Phys. Lett. 
B 579 (2004) 1. 

[8] C.W. Misner, K.S. Thorne and JA. Wheeler, Gravita- 
tion, W.H. Freeman and Co., San Francisco (1973). 

[9] G.L. Alberghi, R. Casadio and A. Gruppuso, Phys. Rev. 
D 61 (2000) 084009. 
[10] R. Brout and G. Venturi, Phys. Rev. D 39, 2436 (1989). 
[11] B.S. DeWitt, Phys. Rev. 160, 1113 (1967); J.A. Wheeler, 
in Batelle rencontres: 1967 lectures in mathematics and 
physics edited by C. DeWitt and J.A. Wheeler, Ben- 
jamin, New York (1968). 



[12] G.L. Alberghi et al, in preparation. 

[13] We choose an operator ordering of convenience since our 
results do not depend on it. For more details, see, e.g., 
R. Steigl and F. Hinterleitner, "Factor ordering in stan- 
dard quantum cosmology," gr-qc/0511149 

[14] A. Tronconi, G. P. Vacca and G. Venturi, Phys. Rev. D 
67, 063517 (2003) 

[15] H.R. Lewis and W.B. Riesenfeld, J. Math. Phys. 10, 1458 
(1969); X.-C. Gao, J.-B. Xu and T.-Z. Qian, Phys. Rev. 
A 44, 7016 (1991). 

[16] F. Finelli, G.P Vacca and G. Venturi, Phys. Rev. D 58, 
103514 (1998). 

[17] N.D. Birrell and P.C.W. Davies, Quantum fields in curved 
space, Cambridge University Press, Cambridge (1982). 

[18] C. Bertoni, F. Finelli and G. Venturi, Class. Quant. Grav. 
13 2375, (1996). 

[19] More general cases will be studied in Ref. [T^. 

[20] For the general theory of invariant operators in quan- 
tum mechanics see Refs. [T^|. and for its application to 
cosmological models see, e.g., Ref. [?ill6|. 



